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Abstract 

We investigate the effect of supersymmetric CP violating phases on the inclusive decay 
B — > Xrf7. Although such a decay contains a large background from B — > Xs^, if isolated 
it may exhibit sizeable CP violation, both in the Standard Model (SM) and in the context 
of models beyond the SM. With unconstrained supersymmetric CP violating phases we 
show that the direct CP asymmetry {Acp) lies in the region —40% < Acp < 40%, where 
a positive asymmetry would constitute a clear signal of physics beyond the SM. Even 
if a direct measurement of i? ^ X^^y proves too difficult experimentally, its asymmetry 
contributes non-negligibly to the measurements of Acp for B Xg^, and thus should be 
included in future analyses. We show that there may be both constructive and destructive 
interference between A'!Pp and A'^^p- 
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1 Introduction 



Theoretical studies of rare decays of b quarks have attracted increasing attention with the recent 
turn-on of the B factories at KEK and SLAC |3];[H- In this paper we are concerned with the 
rare decay b ^ d'j which proceeds via an electromagnetic penguin diagram, and is sensitive to 
the CKM matrix element Vtd- The latter is so far unmeasured directly and in the Wolfenstein 
parameterisation is given by: 

Vtd = AX%l-p-t7]) (1) 

By assuming unitarity of the CKM matrix and using various experimental data {ex, \ Vub/Vcb\, 
Anid, Arris) an allowed region for p and rj at 95% c.l. can be obtained. This enables one to 
infer the range of Vtd consistent with the unitarity of the CKM matrix. 

A direct measurement of Vtd is thus desirable and we will be considering the inclusive decay 
B Xd-f. Current measurements of B^-B^ mixing yield 0.0065 < \VtdV^l\ < 0.010 |], 
where the main source of error is the theoretical uncertainty in the hadronic matrix element 



fsj:^ Bb^- Alternative ways to measure the ratio |Vjd/Vts| have been proposed and include the 
j^oj^o jjiixing ratio Arrid/ Arris and the ratio of the branching ratios (BR) of B ^ Xgl^l^ and 

B^Xdin- §. 

The short distance contribution io B ^ Xd'j is the quark transition b ^7. Estimates 
for the long distance contributions can be found in |P, || and their relative size is expected 
to be around 10% of the short distance contribution. Experimental upper limits exist for the 
branching ratios (BRs) of the exclusive decay channels, B —>■ p^j and B p'^1- CLEO [0 
obtains < 1.7 x 10"^ and < 1.3 x 10~^ respectively, with corresponding measurements by 
BELLE § of < 0.56 x 10"^ and < 2.27 x 10"^ 

There is considerable motivation for calculating the BR and CP asymmetry {Acp) of the 
inclusive channel B Xd'j: 

(i) It provides a theoretically clean way of measuring Vtd, as proposed in ||^, |l^ 



(ii) Acp in the SM is sizeable, and much larger than that for 6 ^ 57 [|10| 



(iii) Acp is sensitive to new physics which contributes to Cj [[IT], |T^, |r^, [T^. 



(iv) The current measurement of Acp for 6 — *• 57 by the CLEO Collaboration |T5[ is sensitive 
to events from b — > dj. Therefore knowledge of Acp for b —>■ dj is essential, in order to 
compare experimental data with the theoretical prediction in a given model. 

We are interested in the effect of unconstrained supersymmetric (SUSY) CP violating phases 
on the inclusive decay BR(5 — >■ Xd'y) . We will be working in the context of the effective SUSY 
model proposed in |]rB|. Such a model allows one to consider the full impact of the phases on the 
rare decays of B mesons, while simultaneously satisfying the stringent bounds on the Electric 
Dipole Moments of the electron and neutron. Such phases may be crucial for generating the 
observed matter-antimatter asymmetry in the universe |]17 . 



Our work is organized as follows. In section 2 we introduce the decays b ^ d'y and b — *■ sj. 
In section 3 we outline our approach to calculate the CP asymmetries, while section 4 presents 
the numerical results. Finally, section 5 contains our conclusions. 



2 The decays b — > d'y and b — > 

Much theoretical study has been devoted to the decay b ^ s^y due to its sensitivity to physics 
beyond the SM [0]. Exclusive channels {B — >■ K*'y etc.) and the inclusive channel have been 



measured at CLEO, ALEPH, BELLE and BaBar |]T9|]. The related decay b ^ d'j has received 
less attention although is expected to be observed at the B factories, at least in some exclusive 
channels. 



Ref. [|l^ calculated BR(5 X^'j) in the context of the SM. It was shown that the ratio R 



defined by 

BR{B ^ X,7) ^ ^ 

is expected to be in the range 0.017 < R < 0.074, corresponding to BR(i? X^'y) of order 
10^^ — > 10^^. With 10^ BB pairs expected from the B factories, one would be able to produce 
10^ —>■ 10^ b d'-f transitions. In the ratio R most of the theoretical uncertainties cancel, and 
hence R may provide a theoretically clean way of extracting the ratio |VJd/^s|- 
The CP asymmetry (Acp), defined hj[] 



is expected to lie in the range —7% < App < —35% in the SM [|T^, where the uncertainty 
arises from varying p and t] in their allowed regions. Also included is the scale dependence (nb) 
of A'qp which occurs from varying mb/2 < /i^ < 2mfe. For definiteness we fix /Zfe = 4.8 GeV, and 
find —5% < App < —28%. Therefore App is much larger than A^p (< 0.6%). By estimating 
values for detection efficiencies, it has been argued in that App may be statistically more 



accessible than A^Jp, at least in the context of the SM. This analysis assumes that B X^'y 
can be clearly isolated from B —>■ Xgj. 

However, it is known that isolating the signal B —>■ X^'j would be an experimental challenge 



since B Xg'y constitutes a serious background. [|I^ has suggested several ways to overcome 
this problem, e.g. demanding a higher energy cut on 7, since 7 from B Xd'y will be more 
energetic than that from B Xg'y. Energy cuts can be used to separate b ^ s'y events from 
charmed background since there is a high photon energy region that is inaccessible to charmed 
states because of the mass of the charm quarks. This method is not feasible for extracting 
6 — i> (i7 events from a 6 — 57 sample. Although the strange quark mass is larger than the 
down quark mass, the respective lightest hadronic single particle final states, K* and p, have 
almost the same mass (they actually overlap strongly). The lightest multi-particle states are 
Kn and tttt, respectively, but even here effects such as bound state effects (neglected in fl^ ) 
smear the spectra out over regions of the order of 200 MeV. These effects constitute one of the 
major theoretical uncertainties in the extraction of BR{B — > ^^7) from the measured part of 
the spectrum, and they make a separation of b d'j and b —>■ s'j via energy cuts impossible. A 
comparison of the photon energy spectra for b ^ sj and b d'-f was made in 0, and showed 
that the photon spectra for both decays are very similar. 



Note that our definition Acp contains a relative minus sign compared to that used in pO|,pl 
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A more promising approach constitute exclusive channels |I|, The improved K/% 

separation at the B factories may enable the inclusive B Xd'~f decay to be reconstructed 
by summing over the relevant exclusive channels as done by CLEO in the measurement of 
B Xg'y. |T^ suggested using a semi-inclusive sample of B ^ ■y + nn decays with a maximum 
of n (say 5) mesons together with a corresponding measurement of-B ^ 'y + K + {n — l)7r. 
The ratio of the widths of the semi-inclusive samples would enable the total inclusive rate to 
be deduced to a very good approximation. Although the extraction of the branching ratio for 
b —>■ dy from exclusive channels might suffer additional uncertainties with respect to b —>■ sj 
the asymmetry should not be affected by these. 

If App and A^Jp cannot be separated, then only their sum can be measured. In the context 
of the SM (with = = 0) the unitarity of the CKM matrix ensures that the sum is zero 



20| , P3[| . This relation holds only for the short distance contribution, which is expected to be 
dominant (c.f. introduction). In the presence of new physics such a cancellation does not occur, 
as will be shown in section 4. As stressed earlier, a reliable prediction of Aqp in a given model 
is necessary since it contributes to the measurement of A^Jp. The CLEO result is sensitive to 
a weighted sum of CP asymmetries, given by: 

Ac^ = + 0.02App (4) 

The latest measurement stands at —27% < ^^tp < (90% c.l.) [0. The small coefficient 
of App is caused by the smaller BR(5 — > Xdj) (assumed to be 1/20 that of BR(5 — > Xsj)) 
and inferior detection efficiencies, but may be partly compensated by the larger value for App. 
We shall see that there can be both constructive and destructive interference between the two 
terms in eq. (|). These effects will be especially important for measurements in future high 
luminosity runs of B factories, in which the precision is expected to reach a magnitude where 
the b ^ dy contribution becomes crucial. For integrated luminosities of 200 fb~^ (2500 fb~^) 



anticipates a precision of 3%(1%) in the measurement of A'^^p 



3 Direct CP Asymmetry in B Xd^sl 

In this paper we explore the effect of CP violating SUSY phases on the direct CP asymmetry of 
the inclusive decay B Xa{s)l- We will show that the asymmetry App may be quite different 
from the SM prediction in a wide region of parameter space consistent with experimental bounds 
from the Electric Dipole Moment (EDM) and BR{B Xsj). 

In our analysis we adapt the "effective SUSY" model, proposed in |]T^ . This model permits 
unrestricted SUSY phases and evades the electric dipole moment (EDM) constraint by invoking 
large masses (of order 20 TeV) for the first two generations of sfermions, thus maintaining their 
contribution to the EDMs within the experimental limits. The third generation sfermions are 
allowed to be relatively light (< 1 TeV) with large phases in /i and the soft breaking term At. 
Such an approach has been used on several occasions |2^, in order to study the maximum 



impact of SUSY phases on rare B decays. It was pointed out in that the third generation 
squarks may contribute to the EDMs via non-negligible two loop diagrams, and we include this 
constraint in our analysis. 

Previous work on the magnitude of App in models beyond the SM include multi-Higgs 
doublet models [jl3|, a Left-Right symmetric model and the MSSM without SUSY 
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phases . A common feature to all these analyses is the possibility of Afjp being of opposite 



sign to that of the SM, which would be a clear signal of new physics. Ref. |]T2[ took At and 



fi to be real, and thus the only source of CP violation was the CKM phase. They found two 
phenomenological acceptable regions, corresponding to 2% < A^p < 21% and -45% < < 
—5%. It is instructive to consider the impact of unconstrained SUSY phases on the inclusive 
decay B Xcfi by taking At and /i complex. The same approach has been used in Ref.[p5 



and it was shown that A^p may lie in the range —16% < A^p < 16%. Other authors p8 
have considered a variety of SUSY models with additional theoretical assumptions, resulting 
in lower values for the maximum value of A^p. The study of b ^ d'y in these models will be 
considered in future work 



We now briefly outline our approach for the calculation of A'^p"' . We assume that flavour 
changing neutral current vertices induced by the gluino and neutralino are absent. Therefore 
to lowest order, the decay at quark level proceeds via the following diagrams, where the photon 
may be emitted from any charged line: 
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t,t 



s, d 



The effective Hamiltonian for b ^ d'y is given by 

AG ^ 

Ueff = — J^v:,Vt,Y.c^{^ib)Q^{^ib) (5) 

where Qi{fib) is the current density operator for the AB = 1 transition and Ci{fib) is its Wilson 
coefficient. The relevant operators for b d'j decay are given by 

Q2 = dL^^CLCLYh, 

Qi = Y^mhdLa^'bRF^,, (6) 



Qs = -^m,dLa^''T%RG%. 



The analogous formulae for the b ^ s'j decay can obtained from Eqs.(|^) and (^) by making 
the replacement d ^ s. 

The asymmetry A'^p"' can be written as: 

A^^P = -—(1.17 X lm[C2C;] - 9.51 X lm[CsC;] + 0.12 x Im[C2C*] 



-9.40 X lm[edis)C2{C; - 0.013C*)] 
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[0.54 Iml^j] - 0.25 Im[^s] - 0.19 Im[C,^s] 



(l + i?e[6])2 + (/m[e7])2 

+3.21 /m[e,(,)(l + 0.65^7* + 0.04^8*)]] (7) 

where ^7,8 = (^7,8 — Cj^)/Cf^ . Cj^s include the total contribution while Cf^ contain only 
the Standard Model one. Here we use Cf*^ = -0.30 and C^^'^ = -0.14. 

In the SM since all the Wilson coefficients are real, the only contribution comes from the 
final term, which corresponds to the CKM phase in ex = V*xVub/V*xVtb- 

The branching ratios in terms of the new physics contributions are given by 

Br(B^X.7) = ^;|^-^|C7rBr(i?^X,ez7,), 

^ (3.48 ± 0.31) X 10-^ [(1 + Re[iy]f + {Im[iy]f] , (8) 

Br(5^X,7) = X'[{l-pf + ri^]Bv{B^X,^). (9) 

The dominant contribution to the decay comes from C7 evaluated at the scale mf,, which 
may be divided into contributions from and respectively: 

Cj = cY + cf + (10) 

Both CY and are purely real, while C7 may possess an imaginary part. We will use the 
leading order expressions for C-j which may be found in . Although higher order corrections 



to CY are available, the corrections to and Cj are only valid in certain limiting cases 
which are not generally applicable [^. Therefore to be consistent to a given order we limit 
ourselves to the leading order expressions for the Wilson coefficients. Such an approach has 
also been adopted in |T^ and [P3|. The magnitude of IC7I is constrained by measurements of 



the branching ratio of 6 ^ 37, and so we only consider points in parameter space that satisfy 
0.2 < |C7(mb)| < 0.38. 

In effective supersymmetry the main contribution to the EDMs of the electron and neutron 
comes from the two-loop constraint induced by Barr-Zee diagrams involving the Higgs pseu- 
doscalar {A^) and the third generation squarks [^. Since the neutron EDM has large hadronic 
uncertainties, we only consider the upper bound of the electron EDM in our analysis. For 
large values of tan (5 this contribution may exceed the present experimental bounds. We have 
checked that the region tan/? < 30 comfortably satisfies this EDM constraint even for a light 
ii (< 200 GeV). In particular, larger values of the Higgs pseudoscalar mass {Ma > 400 GeV) 
induce enough suppression to allow a sizeable parameter space with tan/3 > 30. We shall see 
that the full range of values for .A^p and A^Jp can be obtained for any value of tan jd > 10, and 
so such a constraint has little impact on our results. 

4 Numerical Results 

We vary the (SUSY) parameters in the following range: 
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Here M and M' are respectively the SU{2) and t/(l) gaugino soft masses; At is the tri-hnear 
soft mass for t, with phase 0^; Mq and Mu are the soft masses for the third generation squarks. 
We respect the direct search lower limits on the masses of ti, by discarding generated points 
that do not pass our cuts of m^-^ > 90 GeV and m^± > 80 GeV in addition to the cut on Cj 
mentioned in Sec. 3. We vary p and t] in the range allowed by present CKM fits for the SM 
P]. Note that in the effective SUSY model one should strictly only include the constraint from 
|\4fe/Vcf,|, which corresponds to varying p and r/ in a semi-circular band in the p — 1] plane. 
This enlarged parameter space has little effect on our graphs, except for Fig. 3, which will be 
commented on below. 

If the signal for the inclusive decay can be isolated then a positive asymmetry would be a 
clear sign of new physics. In Fig. |l| we plot App against m^^ , which clearly shows that a light 
ii may drive App positive, reaching maximal values close to +40%. For ii heavier than 250 
GeV the A^p lies within the SM range, which is indicated by the two horizontal lines. 

We note that our upper limit of +40% is larger than the maximum value of 21% attained 
in ||12|. The inclusion of the SUSY phases has joined and expanded the two phenomenological 
regions found in |jl2|, allowing CP asymmetries in the continuous region -40% < A'pp < 40%. 



In Fig. Q we show that the large positive asymmetries can be found anywhere in the interval 



5 < tan/5 < 30, which is the region where the EDM constraint in [27| is comfortably satisfied. 

If the signals from b —>■ sj and b ^ d'y cannot be isolated then one must consider a combined 
signal. In Fig. we plot App against A^p. The maximum values for A^p agree with those 
found in |^ . It can be seen that there is an inaccessible region and the asymmetries can never 
simultaneously be zero e.g. for App ^ 0, \A'Jp\ > 3%. This can be explained from the fact 



that AJp ~ would require Cy to have a sizeable imaginary part in order to cancel the large 
negative contribution from e^. The corresponding effect on A^p would be to cause a sizeable 
deviation from its small SM value. Fig. ^ shows that both A^Jp and App can have either sign, 
resulting in constructive or destructive interference in eq. (|[). If only the \Vub/Vcb\ constraint 
is included in the CKM fits, the enlarged parameter space for p and rj allows much smaller 
asymmetries for App. This is because smaller values of rj are now allowed, which reduces the 
SM contribution to App. The choice of r] ^ would correspond to points in the previously 
inaccessible region. 

In Fig. ^ we plot AF^ + AF^ (defined in eq. (^) against Im(C7). In the SM (as explained 
in Section 2) this sum would be exactly zero in the limit rris = = (neglecting the small 
long distance contribution). From Fig. | it can be seen that AF^ + AF^ is close to if Cy is 
real, the slight deviation being caused by the imaginary parts of the other Wilson coefficients. 
The effect of a non-zero Im(C7) causes sizeable deviations from zero. 

In Fig. H we plot the A^p (defined in eq. (^) against A^p. The right hand plot shows a 
magnification of the area around the origin. The coefficient of App in eq. (^) assumes that 
BR(6 — > d'y)= BR(6 — > s7)/20. Since this ratio of BRs is ~ |Vid/Visp, which in turn is a 
function of the variables p and 77, we replace the factor 1/20 by the above ratio of CKM matrix 
elements. If the contribution from App were ignored in eq. (|^), then Fig. |^ would be a straight 
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line through the origin. The Ajp contribution broadens the hne to a thin band of width ^ 1%, 
an effect which should be detectable at proposed higher luminosity runs of the B factories. 

Note that the width of the line is determined by the amount of b ^ d'y admixture in the 
b s'y sample, eq. (^. In the case of the CLEO measurement the admixture of b ^ d'j 
is about 2.5 times less than the "natural" admixture (ratio of the branching ratios). If the 
experimental analysis can be done with a natural admixture or even ab ^ d'y enriched sample, 
the width of the line would be correspondingly broader. Specifically, for the natural admixture 
the line would be broadened by a factor of 2.5, making b d'-f a 2.5% effect. This effect is the 
same magnitude as the precision attainable with an integrated luminosity of 200 fb~^ at the 
B factories [Q. At this luminosity it will therefore be possible to test the cancellation of the 
asymmetries as predicted by the SM. 

5 Conclusions 

We have studied the effect of supersymmetric (SUSY) CP violating phases on the direct CP 
asymmetry in the inclusive decay B Xd'j. We have performed our calculation in the effective 
SUSY Model, which allows unrestricted SUSY phases without violating the stringent constraints 
from the electron EDM. Although such a decay contains a large background from B — > Xg^, 
it may exhibit large direct CP violation in the context of the Standard Model (SM) and its 
extensions. 

In the SM the CP asymmetry App is expected to lie in the range —5% < Acp < —28%. 
A previous analysis in the context of the MSSM in the absence of SUSY phases found two 
phenomenologically acceptable regions in SUSY parameter space corresponding to asymmetries 
of —45% < Acp < —5% and 2% < Acp < 21%. The latter would constitute a clear signal 
of physics beyond the SM. We have shown that the inclusion of phases in the SUSY breaking 
parameters joins and expands these regions to allow CP asymmetries in the continuous region 
—40% < Acp < 40%, where the exact boundaries depend strongly on the bounds for C7 that 
are imposed. The largest values for the positive asymmetry occur when the stop is lighter than 
200 GeV. Asymmetries of this magnitude are expected to be within the reach of the B factories 
BELLE and BaBar. 

If the inclusive decay B X^'y cannot be isolated from B X^'j then one must consider a 
combined signal. In the SM the sum of the dominant short distance contributions is identically 
zero (for = mg = 0), although such a cancellation does not occur in the effective SUSY 
model since C7 may possess an imaginary part. We have studied the correlation between App 
and Acp, showing that there may be both constructive and destructive interference. The 
contribution of 5 — > Xd'j to the combined signal A'^Jp is suppressed by a branching ratio factor 
l^d/Vtsp, but may be partly compensated by its potentially larger asymmetry. These results 
will be particularly important at proposed high luminosity runs of the B factories, in which 
the experimental error in the measurement of A^Qp is expected to reach a level comparable to 
the magnitude of the B X^'j contribution. 
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